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NUMERICAL  MODEL  OF 
FROZEN  STARTUP  OF 
HIGH-TEMPERATURE  HEAT 
PIPES 

1.1  SUMMARY 

Continuum  and  rarefied  vapor  flows  co-exist  along  the  heat  pipe  length  for  most  of 
the  startup  period.  A  two-region  model  is  proposed  in  which  the  vapor  flow  in  the 
continuum  region  is  modeled  by  the  compressible  Navier-Stokes  equations,  and  the 


I 


vapor  flow  in  the  rarefied  region  is  simulated  by  a  self-diffusion  model.  The  two 
vapor  regions  are  linked  with  appropriate  boundary  conditions,  and  the  heat  pipe 
wall,  wick  and  vapor  flow  are  solved  as  a  conjugate  problem.  The  numerical  solutions 
for  the  entire  heat  pipe  startup  process  from  the  frozen  state  are  compared  with  the 
corresponding  experimental  data  with  good  agreement. 

1.2  INTRODUCTION 

Due  to  the  complex  heat  transfer  phenomena  occurring  within  a  heat  pipe,  numerical 
models  are  normally  used  to  perform  a  detailed  study  of  the  frozen  startup  period. 
Examples  of  previous  works  that  have  developed  numerical  models  include  Chang  and 
Colwell  (1985),  Costello  et  al.  (1987),  Hall  and  Doster  (1990),  Jang  et  al.  (1990),  and 
Cao  and  Faghri  (1993a).  Jang  et  al.  (1990)  described  the  different  startup  periods, 
and  the  general  mathematical  formulation  for  each  period  was  given.  In  an  attempt 
to  solve  the  conservation  equations,  the  effect  of  the  liquid  flow  in  the  porous  wick 
structure  was  neglected,  and  the  vapor  flow  was  simulated  with  a  one- dimensional 
model.  In  the  initial  stage  of  the  simulation,  the  liquid-vapor  interface  was  assumed 
to  be  adiabatic,  and  only  the  heat  conduction  in  the  wall  and  the  melting  in  the  wick 
was  considered.  In  the  transient  continuum  flow  stage,  the  vapor  was  assumed  to  be 
in  a  continuum  flow  state.  An  effort  to  couple  the  continuum  vapor  flow  with  the 
wick  and  wall  was  made,  but  the  coupling  was  apparently  not  completed,  and  the 
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linkage  between  the  two  stages  was  not  considered.  Only  the  model  by  Cao  and  Faghri 
(1993a)  examined  the  entire  heat  pipe  startup  process  from  the  early  free- molecular 
vapor  condition  to  the  final  steady-state  operation.  The  closed-form  analytical  model 
by  Cao  and  Faghri  (1992)  also  determined  the  frozen  startup  characteristics  of  high- 
temperature  heat  pipes.  The  numerical  and  analytical  models  of  Cao  and  Faghri 
(1992,  1993a)  are  described  in  detail  in  this  section. 

Cao  and  Faghri  (1993a)  proposed  a  complete  two-dimensional  numerical 
model  which  split  the  heat  pipe  into  two  regions:  coutinuum  and  rarefied  flows.  In 
the  continuum  flow  region,  the  mass,  momentum,  and  energy  conservation  equations 
were  solved.  In  the  rarefied  flow  region,  the  self-diffusion  model  given  by  Cao  and 
Faghri  (1993b)  was  used.  The  two  regions  were  coupled  with  appropriate  boundary 
conditions  to  model  the  entire  heat  pipe  as  a  conjugate  problem.  This  model  predicted 
the  transient  and  steady  performance  characteristics  of  a  high-temperature  heat  pipe 
started  from  the  frozen  state. 


1.3  RAREFIED  VAPOR  FLOW  REGION 

In  the  early  stages  of  frozen  startup,  the  vapor  pressure  and  density  are  extremely 
low.  Therefore,  the  vapor  is  in  the  rarefied  state,  which  is  somewhat  different  than 
that  in  the  continuum  state.  Also,  a  very  large  vapor  density  gradient  exists  in  the 
axial  direction  of  the  heat  pipe.  Therefore,  the  axial  motion  of  the  vapor  is  caused 
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by  molecular  diffusion.  The  rarefied  vapor  flow  can  be  simulated  by  a  self-diffusion 
model,  where  self-diffusion  refers  to  the  interdiffusion  of  particles  of  the  same  species 
due  to  the  density  gradient.  The  conservation  of  mass  and  energy  equations  for  a 
differential  cylindrical  control  volume  are 
Mass  self- diffusion: 


dp 

dt 


4r(o.¥ 


dz 


dz , 


Energy: 


(1.1) 


i^) + + ;  I  (•*■§ ) + 1  Hr)  o-*) 


where  the  mass  fluxes  Gr  and  G -  are 


Gr  =  pv  =  -Dvfr,  Gz  =  pw  =  -Dv~  (1.3) 

Dv  is  the  coefficient  of  self-diffusion,  which  is  given  by  the  relation  based  on  the 
Chapman-Enskog  kinetic  theory  (Hirschfeider  et  al.,  1966). 


Dv  =  2.628  x  10 


_7 


(1.4) 


Pct2Q,d(kT  /  e) 

where  p  is  the  pressure  in  atmospheres,  a  is  the  collision  diameter  in  A,  e  is  the  max¬ 
imum  energy  of  attraction  between  a  pair  of  molecules,  k  is  the  Boltzmann  constant, 
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M  is  the  molecular  weight,  and  ffc  is  the  collision  integral  for  mass  diffusion.  The 
heat  dissipation  has  been  neglected  due  to  the  low  vapor  temperature  and  small  vapor 
viscosity  during  the  frozen  startup  period. 

1.4  CONTINUUM  VAPOR  FLOW  REGION 

The  vapor  flow  within  heat  pipes  has  been  studied  previously  using  experimental, 
analytical,  and  numerical  methods.  Wageman  and  Guevara  (1960)  simulated  an 
evaporator  by  injecting  air  through  a  porous  pipe  in  order  to  study  the  transition 
to  turbulence.  Their  experimental  results  closely  matched  the  laminar  flow  theory 
by  Yuan  and  Finkelstein  (1956),  even  for  axial  Reynolds  numbers  up  to  10s,  where 
turbulence  should  have  been  well  established.  Bowman  and  Hitchcock  (1989)  studied 
a  simulated  heat  pipe  by  injecting  and  extracting  air  through  a  porous  pipe.  It  was 
found  that  the  vapor  flow  was  always  laminar  in  the  evaporator,  while  it  may  become 
turbulent  in  the  condenser  above  Reynolds  numbers  of  1.2  x  104.  During  the  frozen 
startup  process,  the  axial  vapor  velocity  can  be  supersonic,  but  the  axial  Reynolds 
number  is  relatively  small  due  to  the  low  vapor  densities  of  liquid-metal  working 
fluids.  Faghri  et  al.  (1991b)  showed  that  the  axial  Reynolds  number  was  less  than 
500.  Therefore,  the  continuum  vapor  flow  in  the  heat  pipe  can  be  assumed  to  be 
laminar  along  the  entire  heat  pipe  length. 
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1.5  HEAT  PIPE  WALL  AND  WICK 

The  heat  transfer  through  the  heat  pipe  wall  is  by  conduction.  The  change  of  phase 
of  the  frozen  working  fluid  must  be  considered  in  the  wick  structure,  which  is  as¬ 
sumed  to  be  isotropic  and  homogeneous.  The  phase- change  mode]  by  Cao  and  Faghri 
(1990b)  was  employed  by  Cao  and  Faghri  (1993a,  1993b),  which  is  a  fixed-grid  for¬ 
mulation  having  the  advantage  of  eliminating  time-step  and  grid-size  limitations.  In 
this  analysis,  the  effect  of  the  liquid  flow  in  the  wick  on  the  temperature  distribution 
has  been  neglected.  This  is  due  to  the  fact  that  the  wick  is  usually  very  thin,  the 
thermal  conductivity  of  the  liquid  metal  is  high,  and  heat  conduction  is  the  dominant 
transport  mechanism  in  the  wick  structure.  The  validity  of  this  conduction  model 
in  the  wick  was  also  documented  by  Ivanovskii  et  al.  (1982)  and  Dunn  and  Reay 
(1982).  The  conservation  of  energy  equation  below  is  also  applicable  to  the  entire 
frozen  startup  process. 


l(MCT)  =  l  (t-rg )  -  !<*.« 


(1.5) 


The  coefficients  C  =  C(T),  S  =  S(T),  and  kef[  =  keft(T)  are  given  by 
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C(T)  =  c,ne  +  H/2ST,  (Tmel  -  ST  <T  <  ST  +  Tmcl)  (mushy  phase)  (1-6) 
c/e,  ( T  >  rmei  +  ST)  (liquid  phase) 

Cse(6T  -  T,„el)i  {T  <  ?mel  ~  ST) 

S(T)  =  Cn'ST  +  HI  2  -  (cw  +  H/28T)Tmeh  (TmeJ  -  ST  <T  <  ST  +  Tmd) 

cseST  +  H  —  ceeTlue j,  (T  >  Tme]  4-  <5T) 

(1.7) 


MT)  =  fcse  (fce  -  k,e)  (T  -  T,ne]  +  6T)/28T, 
kie, 


[T  <  rme i  -  *r) 

(rmel  —  st  <t  <  st  +  rmd) 

(r  >  rme,  +  st) 

(1.8) 


where 


Peff  —  ^PPl  4"  (1  lt>)pwst  Cje  —  4”  (1 

Cte  —  4*  ( 1  <p)Cwsi  £ me  —  (Cje  4"  C^p)/2 

For  a  wrapped  screen  wick,  &se  and  A*e  are  calculated  using  the  equation  given  in 


Section  3. 


(1.10) 


k(  [(fcf  4~  k ws)  (I  y)  (k( 

[{ke  +  kws)  +  (1  -  <p)  {kt  -  kw3)} 


1.6  BOUNDARY  CONDITIONS 

A  complete  summary  of  the  boundary  conditions  applicable  to  the  frozen  startup  of 
high-temperature  heat  pipes  can  be  found  in  Table  1.1.  After  the  melting  interface 
reaches  the  wick-vapor  interface,  r  =  /?„,  evaporation  or  condensation  takes  place. 
The  relation  for  the  interface  mass  flux  is 


m"  =  -psi'i 


\  /  M  f  j£_  _  pv  \ 


(1.11) 


where  ms  >  0  for  evaporation,  <  0  for  condensation,  and  a  is  the  accommodation 
coefficient.  p'(  is  the  equilibrium  saturation  pressure  corresponding  to  the  liquid  tem¬ 
perature  Tt,  and  can  be  calculated  using  the  Clausius- Clapeyron  equation. 


Pi 


kfs 

1S\ 

. 

u 

To J. 

(1.12) 


po  and  To  are  the  reference  pressure  and  temperature.  Sukhatme  and  Rohsenow 
(1966)  showed  that  a  may  change  from  0.1  to  1.0,  according  to  the  surface  condition. 

As  was  previously  mentioned,  a  two-region  model  for  the  vapor  flow  in  the 
heat  pipe  was  employed.  The  continuum  and  rarefied  regions  are  coupled  at  the 
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1.1  Summary  of  the  mathematical  models  and  boundary  conditions  for  the 


analysis  of  startup  from  the  frozen  state 


Evaporator  and  condenser  end 
caps,  x  =  0  and  L»: 


Boundary  conditions: 
Heat  pipe  centerline,  r  =  0: 
BZ  -  S£  =0 

Or  ~  Or  V 


Evaporator  and  condenser  end 
caps,  x  =  0  and  Li : 

*-€-o 

Continuum-rarefied  interface^: 

P  =  T  ss  3Vc _ 


t  w,,r  and  ~i,r  are  interfadal  values  from  the  rarefied  vapor  solution. 
XP\,c  and  Ti,c  are  interfacial  values  from  the  continuum  vapor  solution. 


Continuum-rarefied  interface* 
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continuum-rarefied  interface,  which  is  assumed  to  be  planar.  The  two  regions  and  the 
interfacial  grids  are  shown  in  Fig.  1.1. 

It  is  essential  to  maintain  the  mass  balance  across  the  interface  at  each  time 
step.  The  boundary  conditions  at  the  continuum-rarefied  interface  are 

dv 

w  =  Wi,r,  —  =  0,  T  =  71, r  (continuum  vapor  flow)  (1.13) 

UZ 

p  =  /3,jC,  T  =  Ti,c  (rarefied  vapor  flow)  (1-14) 

where  u>,ir  and  Tv  are  the  interfacial  axial  vapor  velocity  and  temperature  on  the 
rarefied  vapor  flow  side,  and  are  obtained  from  the  rarefied  vapor  flow  solution.  p,iC 
and  Ti,r  are  the  interfacial  vapor  density  and  temperature  on  the  continuum  vapor 
flow  side,  and  are  obtained  from  the  continuum  vapor  flow  solution. 

1.7  NUMERICAL  MODELING  AND  PREDIC¬ 
TIONS 

The  conjugate  heat  transfer  and  flow  problem  was  solved  using  the  control-volume, 
finite- difference  technique  with  the  SIMPLE  algorithm.  The  details  of  the  numerical 
scheme  are  given  by  Cao  and  Faghri  (1993a). 
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Continuum  i 
vapor  flow  ■’"j 


Rarefied  vapor  flow 


Continuum-rarefied  interface 


Figure  1.1:  Schematic  of  continuum  and  rarefied  vapor  regions  and  interfacial  grids 


The  frozen  startup  model  was  used  to  simulate  the  early  startup  period  for 
the  experimental  heat  pipe  studied  by  Faghri  et  al.  (1991a,  1991b).  The  multiple- 
evaporator  sodium  stainless-steel  heat  pipe  had  a  vapor  space  radius  R =  10.75  mm, 
wrapped  screen  wick  thickness  6e  =  0.456  mm,  and  pipe  wall  thickness  8W  —  2.15 
mm.  The  details  of  the  heat  pipe  and  experimental  setup  are  given  in  Table  1.2. 
Figure  1.2(a)  shows  the  comparison  between  the  numerical  and  experimental  outer 
wall  temperatures  for  the  early  startup  period  (Case  24a,  Faghri  et  al.,  1991a).  As  can 
be  seen,  the  agreement  is  quite  good  in  the  early  startup  period.  The  accommodation 
coefficient  was  specified  to  be  a  =  0.3,  and  the  phase-change  temperature  range  was 
set  to  2 6T  =  1.0  K  to  simulate  the  phase-change  occurring  at  a  single  temperature. 
The  melting  front  reached  the  wick-vapor  interface  in  the  evaporator  region  quickly 
(about  40  s).  This  means  that  the  startup  period  2  is  very  short. 

The  vapor  Mach  number  for  this  case  is  shown  in  Fig.  1.2(6).  Since  the 
active  evaporator  was  located  away  from  the  evaporator  end  cap,  the  vapor  flows  both 
toward  the  condenser  and  the  evaporator  end  caps.  Very  early  in  the  startup,  the 
vapor  density  was  very  low,  and  the  Mach  number  was  high.  As  vapor  accumulated 
in  the  active  heater  section  due  to  evaporation,  the  vapor  velocity  decreased  to  the 
subsonic  state  (Ma  <  1 )  due  to  the  larger  vapor  density.  Away  from  the  active  heater, 
the  vapor  density  was  very  low,  and  the  vapor  velocity  reached  the  supersonic  state 
in  both  directions.  The  axial  vapor  density  distribution  is  shown  in  Fig.  1.2(c). 
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Table  1.2:  Design  summary  of  the  high-temperature  heat  pipe  with  multiple  heat 


sources 


Container  material 

- -  ~  ■ 

304L  Stainless  Steel 

Wick  material 

304  Stainless  Steel 

End  cap  material 

304L  Stainless  Steel 

Length 

1000  mm 

Container  o.d. 

26.7  mm 

Container  i.d. 

22.4  mm 

End  cap  thickness 

3.175  mm 

Screen  mesh  number 

100  mesh  (3.94  mm-1) 

Screen  wire  diameter 

0.114  mm 

Screen  wick  thickness 

0.456  mm 

Vapor  core  diameter 

21.5  mm 

Working  fluid 

sodium 

Fluid  charge 

30,  43  g 

Evaporator  length* 

(4)  @  53.0  mm  each 

Transport  length 

188  mm 

Condenser  length** 

292  mm 

(*)=Evaporators  separated  by  90  mm  adiabatic  sections  with  a  20  mm 

adiabatic  section  between  the  evaporator  end  cap  and  evaporator  1 

(**)=20  mm  adiabatic  section  between  the  condenser  and  the  condenser  end  cap 


Figure  1.2:  Heat  pipe  frozen  startup  in  the  early  stages  (Case  24a,  evaporator  2  active, 
Q  =  60  W):  (a)  Outer  wall  temperature  compared  with  experimental  data  of  Faghri 
et  al.  (1991a);  (6)  Vapor  Mach  number;  (c)  Vapor  density 
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Although  the  density  in  the  active  heater  region  increased  significantly  with  time, 
very  large  density  gradients  existed  in  both  directions,  and  the  density  was  very  low 
outside  the  active  evaporator  region.  Solving  eqn.  ()  for  the  vapor  density  when 
continuum  flow  is  established  results  in  pv  >  3.3  x  10“4  kg/m3.  From  Fig.  1.2(c), 
it  can  be  seen  that  continuum  flow  is  reached  after  810  s  only  in  the  active  heater 
section,  while  the  vapor  flow  in  the  rest  of  the  heat  pipe  is  in  the  rarefied  state. 

As  discussed  in  Section  2,  a  may  change  significantly  depending  on  the 
surface  condition,  and  there  is  little  information  available  concerning  its  value  for  the 
wick-vapor  interface  of  a  heat  pipe.  The  influence  of  a  on  the  numerical  simulation 
was  examined  by  running  Case  24a  again  with  different  values  of  a,  as  shown  in 
Fig.  1.3.  As  can  be  seen,  the  numerical  results  for  the  outer  wall  temperature  are 
insensitive  to  a.  Since  the  wall  temperature  depends  on  the  interfacial  mass  flux  m'g', 
which  in  turn  depends  on  both  a  and  the  pressure  difference  between  the  wick-vapor 
interface  and  the  vapor  flow,  a  larger  value  of  a  resulted  in  a  higher  vapor  pressure  and 
reduced  the  pressure  difference  between  the  wick-vapor  interface  and  the  vapor  flow. 
As  a  result,  the  two  factors  offset  each  other,  so  the  interfacial  mass  flux  remained 
relatively  constant. 

Figure  1.4  compares  the.  numerical  and  experimental  outer  wall  tempera¬ 
tures  during  all  of  the  frozen  startup  periods  for  Cases  1  la-1  If  (Faghri  et  al.,  1991a) 
with  evaporator  1  active,  which  is  closest  to  the  evaporator  end  cap.  The  agreement 
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Figure  1.3:  Outer  wall  temperatures  with  different  accommodation  coefficients  at  t 
=  675  s  for  Case  24a 


Figure  1.4:  Outer  wall  temperatures  compared  with  experimental  data  for  Cases 
lla-llf 


17 


between  the  numerical  solutions  and  the  experimental  data  is  generally  good.  The 
total  heat  inputs  in  the  different  cases  were  maintained  at  different  levels  in  the  ac¬ 
tive  evaporator  and  incremented  from  311  W  (Case  11a)  to  598  W  (Case  Ilf).  The 
corresponding  power  outputs  at  the  condenser  were  119  W  (Case  11a)  and  298  W 
(Case  Ilf).  The  hot  zone  moved  from  the  active  evaporator  down  the  length  of  the 
heat  pipe  with  time.  At  t  =  185  min,  which  corresponds  to  the  steady  state  for  Case 
11a,  the  hot  zone  reached  the  condenser  section,  but  the  low  heat  input  prevented  the 
hot  zone  from  reaching  the  condenser  end  cap.  At  t.  >  185  min,  the  heat  input  was 
increased  to  the  higher  level  in  Case  lib,  and  the  wall  temperature  increased.  At  the 
heat  input  corresponding  to  Case  11c,  the  hot  zone  filled  the  heat  pipe  completely, 
and  continuum  vapor  flow  was  established  along  the  entire  heat  pipe  length.  Cases 
lla-llc  covered  the  frozen  startup  process.  After  Case  11c,  the  performance  is  the 
transient  heat  pipe  operation  without  phase  change  from  the  frozen  state.  For  the 
numerical  results  at  the  steady  state,  the  total  heat  output  at  the  condenser  was 
within  1.5%  of  the  input  at  the  evaporator  in  all  cases. 

The  double-wall  artery  sodium  heat  pipe  studied  by  Ponnappan  (1989)  was 
then  simulated.  The  heat  pipe  was  2  m  in  length,  with  a  vapor  radius  of  Rv  =  6.35 
mm,  an  outer  wall  radius  of  R0  =  11.1  mm,  and  a  stainless  steel  wall  of  thickness 
Sw  =  1.65  mm.  The  details  of  the  experimental  setup  are  given  in  Table  1.3.  Fig¬ 
ure  1.5  shows  the  comparison  between  the  numerical  and  experimental  outer  wall 
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Table  1.3:  Design  summary  of  the  gas- loaded  high-temperature  double-wall  artery 
heat  pipe  (Ponnappan,  1989) 


Outer  tube,  inner  tube,  end  caps 

304  Stainless  steel 

Screen  wick  material 

304  Stainless  steel 

Fill  valve 

Stainless  steel  bellows  valve 

with  stellite  tip  stem 

Working  fluid 

sodium 

Noncondensible  gas 

argon 

Temperature  range 

400  to  850°  C 

Initial  temperature 

25°  C 

Heat  transport  capacity 

1800  W  at  800°  C 

Effective  length 

1.388  m 

Total  length 

2.03  m 

Evaporator/adiabatic/  condenser 

0.375/0.745/0.91  xn 

Overall  diameter 

2.22  x  10'2  m 

Vapor  core  diameter 

1.27  x  lO"2  m 

Fluid  charge 

92.04  g 
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Figure  1 .5:  Outer  wall  temperatures  compared  with  experimental  data  given  by  Pon- 
nappan  (1989)  for  Q  =  289.6  W 


temperatures  for  a  heat  rate  of  Q  —  289.6  W.  Again,  the  agreement  is  generally  good. 

Figure  1.6  shows  a  comparison  between  the  numerical  vapor  and  wall  tem¬ 
peratures  for  Case  11a.  For  most  of  the  previous  frozen  startup  studies,  the  wall 
temperatures  were  obtained  from  experimental  measurements,  and  the  vapor  temper¬ 
atures  were  assumed  to  be  of  the  same  magnitude  and  trend  as  the  wall  temperatures. 
This  is  not  entirely  correct,  however,  especially  in  the  slope  region,  as  shown  in  Fig. 
1.6.  The  vapor  temperature  front  is  always  steeper  and  is  located  close  to  the  middle 
of  the  slope  region  of  the  wall  temperature  distribution.  In  order  for  the  hot  zone 
to  move  forward,  heat  must  be  supplied  by  the  vapor  to  the  slope  region  to  increase 
the  wall  and  wick  temperature  by  condensation  in  that  region.  Therefore,  a  large 
temperature  difference  between  the  wall  and  vapor  is  required. 

Figure  1.7  shows  the  axial  Mach  number  along  the  heat  pipe  for  Cases  11a- 
1  If.  The  Mach  number  increased  in  the  evaporator,  decreased  close  to  the  continuum- 
rarefied  interface,  and  then  increased  again  in  the  rarefied  region  before  dropping  to 
zero  near  the  condenser  end  cap.  In  the  evaporator,  the  Mach  number  began  to  rise 
due  to  evaporation  at  the  wick-vapor  interface.  After  the  exit  of  the  evaporator,  the 
axial  vapor  mass  flow  rate  decreased  due  to  condensation.  However,  the  vapor  density 
also  decreased  significantly  in  this  region.  The  rate  at  which  the  density  decreased 
exceeded  the  decrease  of  the  axial  mass  flux  due  to  condensation.  Therefore,  the 
vapor  velocity  and  Mach  number  increased  further.  When  the  vapor  approached  the 


21 


Figure  1.6:  Comparison  of  centerline  vapor  and  wall  temperatures  for  Case  1  la 
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hot  zone  front  of  the  wall  temperature,  intense  condensation  occurred.  The  Mach 
number  decreased  to  a  much  lower  value  before  entering  the  rarefied  region.  The 
very  large  Mach  number  in  the  rarefied  flow  is  solely  due  to  the  extremely  low  vapor 
density.  Although  the  vapor  velocity  is  very  high,  the  mass  flux  is  very  low. 
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CLOSED-FORM  ANALYTICAL 
MODEL  OF  FROZEN  STARTUP 
OF  HIGH-TEMPERATURE 
HEAT  PIPES 


2.1  SUMMARY 


The  previous  numerical  and  experimental  studies  indicate  that  the  high  tempera¬ 
ture  heat  pipe  startup  process  is  characterized  by  a  moving  hot  zone  with  relatively 
sharp  fronts.  Based  on  the  above  observation,  a  fiat-front  model  for  an  approximate 
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analytical  solution  is  proposed.  A  closed-form  solution  related  to  the  temperature  dis¬ 
tribution  in  the  hot  zone  and  the  hot  zone  length  as  a  function  of  time  are  obtained. 
The  analytical  results  agree  well  with  the  corresponding  experimental  data,  and  pro¬ 
vide  a  quick  prediction  method  for  the  heat  pipe  startup  performance.  Finally,  a 
heat  pipe  limitation  related  to  the  frozen  startup  process  is  identified,  and  an  explicit 
criterion  for  the  high  temperature  heat  pipe  startup  is  derived.  The  frozen  startup 
limit  identified  in  this  report  provides  a  fundamental  guidance  for  high  temperature 
heat  pipe  design. 

2.2  INTRODUCTION 

In  the  previous  section,  the.  frozen  starts  period  of  high-temperature  heat  pipes 
was  successfully  simulated  by  a  comprehensive  detailed  numerical  model.  The  model 
described  the  entire  frozen  startup  process,  and  compared  well  with  experimental 
data.  While  the  numerical  model  is  robust  and  is  able  to  solve  the  problem  in  a  general 
way,  the  time  needed  to  develop  the  model  and  the  computer  resources  needed  to  run 
the  program  make  the  model  costly.  Therefore,  an  analytical  model  is  presented  here, 
so  that  important  system  parameters  can  be  related.  While  the  applicability  of  such 
an  analytical  solution  may  be  limited,  it  can  aid  in  the  understanding  of  the  physical 
phenomena  occurring  in  the  problem. 
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Cotter  (1967)  developed  an  analytical  model  which  incorporated  the  frontal 
startup  behavior  with  a  lumped  one-dimensional  effective  thermal  conductivity  model. 
A  major  limitation  of  the  above  model  was  the  need  to  know  what  the  temperature 
of  the  heat  pipe  will  be  at  a  given  power  input,  or  what  power  will  be  transported  at 
a  given  temperature  before  using  the  analytical  results.  For  heat  pipe  operation,  the 
relation  between  the  working  temperature  and  the  power  input  is  unknown  a  priori. 

Sockol  and  Forman  (1970)  observed  that  at  sufficiently  high  heat  input  lev¬ 
els,  the  evaporator  wall  temperature  rose  to  some  intermediate  level  during  frozen 
startup  of  a  lithium  heat  pipe,  and  remained  almost  constant  as  a  steep  temperature 
front  moved  down  the  length  of  the  pipe.  When  the  hot  zone  reached  the  con¬ 
denser  end  cap,  the  heat  pipe  temperature  increased  to  a  steady-state  value.  Based 
on  these  observations,  Sockol  and  Forman  assumed  the  hot  zone  wall  temperature 
was  uniform,  and  the  cold  zone  remained  at  the  ambient  temperature.  Two  lumped 
first-order  differential  equations  for  the  wall  temperatures  were  derived  by  making  an 
energy  balance  over  the  evaporator  and  condenser  sections.  A  one-dimensional  vapor 
flow  model  was  used,  and  it  was  assumed  that  the  vapor  velocity  was  sonic  at  the 
end  of  the  hot  zone.  The  temperature  of  the.  hot  zone  was  obtained  by  solving  the 
four  coupled  equations  numerically.  The  numerical  hot  zone  temperatures  were  com¬ 


pared  with  experimental  data,  and  were  50  to  80° C  higher  than  the  measured  values. 
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Since  the  uncertainty  of  the  measured  data  was  ±30° C,  the  numerical  solution  was 
considered  to  be  in  good  agreement. 

Cao  and  Faghri  (1992)  derived  a  closed-form  analytical  solution  based  on 
the  flat-front  startup  model.  The  results  of  the  analytical  solution  are  compared  in 
the  following  section  to  experimental  data  and  the  numerical  solutions  given  by  Cao 
and  Faghri  (1993a)  with  good  agreement. 

2.3  ANALYTICAL  SOLUTION 

At  sufficiently  high  heat  input  levels,  the  evaporator  wall  temperature  of  a  liquid- 
metal  heat  pipe  will  rise  to  some  intermediate  level,  and  remain  almost  constant  as 
a  steep  temperature,  front  moves  down  the  length  of  the  pipe  (Sockol  and  Forman, 
1970;  Faghri  et  ah,  1991a;  Ponnappan,  1989).  The  numerical  solutions  in  the  previ¬ 
ous  section  showed  that  the  vapor  hot  zone  front  was  steeper  than  that  of  the  wall 
temperature,  and  was  located  nearly  in  the  middle  of  the  slope  region  of  the  wall 
temperature  distribution.  Based  on  these  facts,  a  flat-front  frozen  startup  model  is 
proposed  to  analyze  the  startup  process  of  high-temperature  heat  pipes.  At  time  t 
and  t  4-  St,  the  actual  heat  pipe  wall  temperature  distributions  are  shown  schemati¬ 
cally  in  Fig.  2.1  by  two  solid  curves.  The  wall  temperatures  were  approximated  by 
the  two  dashed  curves  with  linear  profiles,  where  flat  fronts  are  located  in  the  middle 
of  the  slope  regions  of  the  wall  temperature  distributions.  The  temperatures  of  the 
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flat  fronts  were  assumed  to  be  at  the  continuum-rarefied  transition  temperature,  jTtr. 
The  hot  zone  wall  temperature  at  the  evaporator  end  cap  is  7)2.  The  hot  zone  lengths 
at  time  t  and  t  +  St  are  l  and  1  +  81 ,  respectively,  where  both  T/t  and  l  are  functions 
of  time. 

The  temperature  of  the  evaporator  section  (0  <  z  <  Lt)  increases  as  heat  is 
applied.  Therefore,  the  hot  zone  length  1  should  be  equal  to  Le  before  the  temperature 
at  the  exit  of  the  evaporator  reaches  the  transition  temperature  Ttr  and  the  hot  zone 
begins  to  move  into  the  adiabatic  section.  For  l  —  Le,  a  lumped  energy  balance  results 
in  an  expression  for  the  time  needed  for  the  evaporator  exit  temperature  to  reach  the 
transition  temperature  when  Q  is  constant. 


t.  =  CLt(Taw  -  IQ  (2.1) 

where  C  is  the  heat  capacity  per  unit  length  of  the  heat  pipe,  which  includes  the  heat 
capacities  of  the  wall  and  wick,  and  T is  the  initial  (ambient)  temperature  of  the 
heat  pipe.  The  average  hot  zone  wall  temperature,  Taw,  is  based  on  a  linear  profile 

Taw  =  (Tk  +  Ttr)/ 2  =  AT/2  +  Ttr  (2.2) 

where  AT  =  Th  —  Tlr  is  the  total  temperature  drop  across  the  hot  zone. 

An  energy  balance  over  the  time  interval  t  and  t  +  6t  in  terms  of  average 
hot  zone  temperature  Taw  for  the  hot  zone  wall  and  wick  gives 
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(2.3) 


C~  [l(Taw  -  T0 o)]  +  h"  [l(Tau,  -  T^)]  =  Qm 

where 

h*  =  2tt R0h,  Qm  =  Q  +  -  Too)(Le  +  L0)  (2.4) 

Q  is  the  evaporator  heat  input,  and  h  is  the  heat  transfer  coefficient  at  the  condenser. 
The  radial  temperature  difference  in  the  wick  and  the  wall,  and  the  heat  capacity 
of  the  vapor  have  been  neglected  in  the  above  formulation.  For  noncircular  heat 
pipes,  2nR0  should  be  replaced  by  the  circumferential  length  W.  In  many  cases, 
frozen  startup  is  accompanied  by  a  radiative  boundary  condition  at  the  condenser. 
Therefore, 


h'  =  Trrr;  =  +  r„)  (2.5) 

where  qc  is  the  heat  flux  from  the  surface  of  the  condenser.  Equation  (2.3)  has  a 
general  solution  assuming  h*  is  constant. 


l{Taw  -  Too)  =  exp  [  (^r)  exp  dt  +  cx  exp  ^  (2.6) 

where  U  is  the  initial  time  from  which  the  solution  is  to  be  considered,  and  cx  is  a 
constant  to  be  determined  by  the  initial  condition.  In  the  adiabatic  section,  Le  <  l  < 
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[Lt  +  La)(T{ 


It  has  been  shown  experimentally  that  the  average  hot  zone  temperature  Taw  is  a 


weak  function  of  time  for  a  constant  Q.  Therefore,  Q *  can  be  assumed  to  be  constant 
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in  the  integration,  and  eqn.  (2.10)  can  be  simplified  to 


l  = 


Q 


+  {Le  +  La) 


[hm(Taw  -  r^) 

( Lf  -f  Lg)(Tgw<i  —  Tpg) 

{TgW  -  Too) 


1  -  exp  — 


exp 


h'(t  -  tj)\ 

C  ) 


C 


(2.11) 


Equation  2.11  is  applicable  for  /  >  (Lf.  +  La),  and  eqns.  (2.1)  and  (2.7)  should  be 
used  for  l  =  Le  and  Le  <  l  <  (Le  +  La),  respectively.  It  can  be  seen  that  the 
above  equations  are  functions  of  time  t  and  average  hot  zone  wall  temperature  Taw. 
Therefore,  the  hot  zone  temperature  distribution  is  found  in  the  next  section. 

The  heat  pipe  wall  temperature  is  directly  related  to  the  vapor  temperature 
distribution.  By  neglecting  the  radial  heat  conduction  in  the  wick  and  wall,  the 
temperature  drop  of  the  wall  in  the  hot  zone,  AT,  can  be  approximated  to  be  the 
same  as  that  of  the  vapor.  It  was  shown  that  the  total  pressure  drop  for  the  laminar 
flow  of  vapor  along  the  length  of  a  heat  pipe  is 


=  rrfnr  ^ ^ 

K  pxillyh'fg 

For  heat  pipes  without  an  adiabatic  section,  La  =  0,  eqn.  (2.12)  reduces  to 


(2.12) 


AP=-^-(i.  +  «=  4"Qi’ 


(2.13) 


*P«Rlhh  ”P,.Rthh 

where  Lt  is  the  total  heat  pipe  length.  In  the  flat-front  startup  model,  the  flat  front 
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at  the  interface  between  the  hot  and  cold  zones  is  equivalent  to  the  condenser  end 
cap.  Therefore,  eqn.  (2.13)  is  applied  only  to  the  vapor  flow  in  the  hot  zone  with  l 
replacing  Lt. 


A  p  = 


4/ iQl 


(2.14) 


*PvRthf9 

The  temperature  drop  AT  corresponding  to  the  pressure  drop  A p  is  found  by  inte¬ 
grating  the  differential  form  of  the  Clausius-Clapeyron  equation 


AT  = 


A  p 


(2.15) 


where  (T/  pvhjg)av  is  the  average  value  over  the  integration  interval.  The  temperature 
drop  can  be  found  by  combining  eqns.  (2.14)  and  (2.15). 


AT  = 


4  TavpQl 


(2.16) 


*PlRlh% 

To  account  for  vapor  compressibility  and  density  variation,  the  vapor  density  in  eqn. 
(2.16)  is  assumed  to  be  the  average  density  corresponding  to  Tav.  Combining  the 
ideal  gas  law  and  the  integrated  form  of  the  Clausius-Clapeyron  equation  gives 


Po 

Pv  =  p  7p-  eXP 

ai' 


JHi  ( _L  _  IS] 

Ra  \Tav  To ) 


(2.17) 


where  po  and  To  are  the  reference  (saturation)  pressure  and  temperature.  By  assuming 
a  linear  hot  zone  temperature  and  neglecting  the  radial  temperature  drop  across  the 
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wall  and  wick,  we  find  that  Tav  =  Taw,  where  Taui  is  given  by  eqn.  (2.2).  Combining 
eqns.  (2.16),  (2.17)  and  (2.2)  gives 

h,,(  1  1  \  \n,(T,r  +  A372)  (MT^  +  AT/Z^QlV12] 

R,  \T„  T„  +  AT/2J  [  p„  V  )  j 

(2.18) 

Equation  2.18  can  be  solved  iteratively  using  the  Newton-Raphson/secant  method 
for  the  temperature  drop,  AT,  and  Taw  and  Tk  can  be  obtained  from  eqn.  (2.2). 

The  analytical  closed-form  solutions  for  the  hot  zone  length  l  as  a  function 
of  time  for  the  frozen  startup  period  are  obtained  from  eqn.  (2.1)  for  l  =  Le,  eqn. 
(2.7)  for  Le  <  l  <  (Le  +  Ta),  and  eqn.  (2.11)  for  (Le  +  La )  <  l  <  Lt.  In  all  cases, 
eqn.  (2.18)  is  used  to  find  AT,  and  eqn.  (2.2)  is  used  to  determine  Tau,  and  T^.  The 
equations  are  solved  sequentially  as  follows:  For  an  initial  hot  zone  length  of  /,  at  an 
initial  time  t,-,  the  temperature  drop  over  the  hot  zone  is  found  from  eqn.  (2.18).  The 
average  hot  zone  wall  temperature  Taw  and  the  evaporator  end  cap  temperature  Th 
are  obtained  from  eqn.  (2.2).  Then,  the  hot  zone  length  at  time  i,-  -f  AT  where  At 
is  the  time  step,  can  be  calculated  using  eqn.  (2.7)  or  eqn.  (2.11),  depending  on  the 
location  of  the  hot  zone  front.  The  values  of  AT,  Taw,  and  Th  are  calculated  at  time 
t,  -f  At  using  the  corresponding  hot  zone  length  obtained. 

The  closed-form  analytical  solution  was  compared  to  the  experimental  data 
for  Case  13  of  Faghri  et  al.  (1991a,  1991b)  in  Fig.  2.2.  The  solution  started  from 
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Figure  2.2:  Analytical  wall  temperatures  compared  with  experimental  data  for  Case 


the  inlet  of  the  condenser,  and  eqn.  (2.11)  was  used  to  calculate  l  from  TaU).  The 
agreement  between  the  analytical  results  and  the  experimental  data  is  generally  good. 
Figure  2.3  also  compares  the  analysis  with  data  from  the  case  with  Q  =  289.6  W  from 
Ponnappan  (1989).  The  analytical  results  agree  very  well  with  the  experimental  data, 
except  the  curves  nearest  the  condenser  end  caps.  The  data  in  that  location  indicated 
that  vapor  flow  was  still  in  the  rarefied  state.  In  eqn.  (2.16),  it  can  be  seen  that  the 
temperature  drop  over  the  hot  zone  is  dependent  on  the  fourth  power  of  the  vapor 
space  radius,  Rv.  For  the  heat  pipe  of  Faghri  et  al.  (1991a),  R„  =  10.75  mm,  and 
the  temperature  drop  over  the  hot  zone  was  relatively  small.  For  the  heat  pipe  by 
Ponnappan,  Rv  —  6.35  mm,  and  the  temperature  drop  was  significant.  The  analytical 
closed-form  solution  predicted  the  temperature  drops  over  the  hot  zones  for  the  two 
heat  pipes  remarkably  well. 

Figure  2.4  compares  the  numerical  and  analytical  solutions  of  the  outer  wall 
temperature  with  the  experimental  data  for  Case  1  la  of  Faghri  et  al.  (1991a).  The  flat 
fronts  of  the  analytical  results  intercept  the  numerical  curves  almost  in  the  middle  of 
the  slope  regions,  and  compare  well  with  the  numerical  solution  and  the  experimental 
data. 
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Figure  2.3:  Analytical  wall  temperatures  compared  with  experimental  data  given  by 
Ponnappan  (1989)  for  Q  =  289.6  W 


38 


39 


2.4  FROZEN-STARTUP  LIMITATION 

Heat  pipe  designs  are  generally  based  on  the  well-known  performance  limitations 
which  are  discussed  extensively  in  Section  4.  These  are  the  sonic  vapor  velocity  limit, 
the  wick  boiling  limit,  and  the  wick  capillary  limit.  During  the  frozen  startup  process, 
the  heat  pipe  may  encounter  the  sonic  limit,  but  this  is  a  temporary  condition  which 
will  disappear  when  the  heat  pipe  working  temperature  rises  to  a  sufficiently  high 
level.  The  boiling  limit  rarely  occurs  in  high-temperature  heat  pipes  due  to  the  high 
thermal  conductivity  of  the  liquid-metal  working  fluid  and  the  thin  wick  structure. 
The  capillary  limit  usually  will  not  be  encountered  during  the  frozen  startup  process 
since  the  active  length  of  the  heat  pipe  is  less  than  the  total  length.  Therefore,  with 
a  proper  design,  the  capillary  limit  should  not  occur.  However,  when  heat  is  applied 
only  in  the  evaporator  region  of  a  frozen  heat  pipe  during  startup,  the  working  fluid 
may  be  frozen  out  of  the  evaporator  to  the  adiabatic  or  condenser  regions.  In  this 
case,  a  frozen- start  up  limit,  which  is  different  than  the  other  limitations,  needs  to  be 
considered  (Cao  and  Faghri,  1992). 

The  frozen  startup  process  is  shown  schematically  in  Fig.  2.5.  The  input 
heat  flux  qe  is  applied  at  the  evaporator  and  evaporation  takes  place  at  the  wick-vapor 
interface.  The  vapor  flows  down  the  active  length  and  condenses  on  the  wick  surface, 
releasing  the  latent  heat  of  vaporization.  A  liquid-solid  interface  in  the  wick  moves 
down  the  heat  pipe  with  velocity  Wi  due  to  condensation  onto  the  wick.  The  vapor 
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liquid 

state 


leaves  the  evaporator  section  with  mass  flow  rate  me.  is  part  of  the  total  mass 
flow  rate  that  condenses  on  the  wick  where  the  working  fluid  is  melted.  This  part 
supplies  liquid  condensate  to  be  pumped  back  to  the  evaporator  by  the  capillary  wick. 
The  part  of  the  mass  flow  rate  of  vapor  which  condenses  on  the  solid-wick  region,  mcl, 
is  frozen  out  and  cannot  flow  back  to  the  evaporator  section.  However,  the  movement 
of  the  liquid-solid  interface  in  the  wick  down  the  length  of  the  heat  pipe  will  increase 
the  mass  of  the  liquid  working  fluid,  which  compensates  for  the  mass  loss  due  to  the 
freezing  out  of  vapor  onto  the  frozen  solid-wick  surface.  A  mass  balance  of  the  liquid 
working  fluid  in  the  liquid-wick  region  requires 

< PPtA'wWi  >  m(,  -  met  =  mcs  (2.19) 

where  <p  is  the  porosity  of  the  wick,  pt  is  the  density  of  the  liquid,  and  A'w  is  the  cross- 
sectional  area  of  the  wick  saturated  by  the  working  fluid.  Equation  (2.19)  indicates 
that  the  increase  in  the  mass  of  liquid  in  the  wick  due  to  the  liquid-solid  interface 
movement  must  be  greater  than  or  equal  to  the  mass  loss  due  to  the  vapor  freezing 
onto  the  solid-wick  interface.  If  this  condition  is  not  met,  liquid  depletion  in  the 
evaporator  occurs  due  to  the  vapor  being  frozen  out.  This  relation  was  observed  and 
defined  as  the  frozen-startup  limit  (FSL)  originally  by  Cao  and  Faghri  (1992). 

Frozen  startup  is  characterized  by  the  movement  of  the  hot  zone  down  the 
length  of  the  heat  pipe.  At  time  t ,  as  shown  in  Fig.  2.6,  the  location  oJ  the  melting 
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Figure  2.6:  Description  of  the  criterion  for  the  frozen-startup  limit 
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front  is  at  z  =  l.  At  t  +  8t,  the  melting  front  moves  to  1  +  81.  The  melting  temperature 
Tmei  intercepts  the  two  temperature  distributions  at  points  a  and  d.  The  heat  energy 
received  by  the  part  of  the  slope  region  whose  temperature  is  below  Tmei  during  the 
time  interval  8t  is 


Qc,8t  =  C  ■  Aa-b-c.d  (2.20) 

where  C  is  the  heat  capacity  per  unit  length  and  An-b-c-d  is  the  area  enclosed  by  the 
curve  connecting  points  a,  6,  c,  and  d.  Since 

~  v4a_(,-_C'_rf  =  (:Tnici  —  Too)8l  (2-21) 

eqn.  (2.20)  can  be  approximated  to  be 

QcJt  =  C(Tmd  -  T«o)«  (2.22) 

The  corresponding  mass  of  vapor  that  condensed  onto  the  frozen  wick  during  8t, 
neglecting  axial  conduction,  is 

8mcs  =  mcJt  =  QcJt/hfg  =  C(Tmci  -  Too)8l/hfg  (2.23) 

The  mass  condensation  rate  is 
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c  /fj  C(Tlnd-  Too)  61 

rncs  =  6mcs/6i  = - - - —  (2.24) 

kjg  6t 

The  rate  of  increase  of  the  liquid  mass  is  ^ppi A'JW^  where  the  time  rate  of  change  of 
the  hot  zone  length  is 

ft  =  *  <2-25> 

Substituting  eqns.  (2.24)  and  (2.25)  into  eqn.  (2.19)  and  rearranging  gives 

FSL  =  >  1  (2.26) 

Equation  (2.26)  is  an  approximate  expression  of  the  frozen-startup  limit  for  high- 
temperature  heat  pipes.  For  a  circular  cross-section  heat  pipe, 

A'w  =  2irRv8'e  (2.27) 

where  R \  is  the  vapor  space  radius  and  8\  is  the  thickness  of  the  working  fluid  in  the 
wick  structure.  It  should  be  pointed  out  that  for  normal  operation,  8't  may  not  be 
equal  to  the  wick  structure  thickness  8 f.  Also,  for  a  circular  heat  pipe, 

C(Tmti  -  T^)  ~  (2 xRySfCintPofi  +  2*R06wcwPw){Tmi  -  T^)  +  2x Rv8'eHpeip  (2.28) 
where  Cm,e  is  the  average  effective  specific  heat  of  the  wick  over  the  liquid-wick  and 
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solid-wick,  pen  is  the  effective  density  of  the  wick,  8W  is  the  wall  thickness,  and  H  is 
the  latent  heat  due  to  melting.  Equation  (2.26)  can  he  rewritten  as 


FSL  = 


(pP(RVSfhfg 


>  1 


(2.29) 


(ftv^iCrnePeft  “h  R-o^ui^ui Pui ) ( -Fmel  Too)  "f"  Ry8(H Pl^p 

For  a  thick  wall  or  wick,  R0  and  R ^  in  eqn.  (2.29)  should  be  replaced  by  the  average 
radius  of  the  wall,  Rwo  =  (Rw  -1-  R0)/ 2,  or  the  average  radius  of  the  wick,  ^  = 


( Rv  +  R%v)/ 2. 

The  frozen-startup  limit  v.  as  applied  to  the  heat  pipe  studied  by  Faghri  et 
al.  (1991a),  which  resulted  in 


FSL  ~  1.46  >  1 

This  indicates  that  the  heat  pipe  should  not  have  encountered  the  frozen-startup 
limit,  which  was  indeed  the  case.  It  should  be  noted  that  FSL  was  very  close  to  unity 
for  the  heat  pipe  studied  by  Faghri  et  al.,  which  means  that  the  frozen-startup  limit 
does  need  to  be  considered  in  the  design  process.  The  calculation  was  also  made  for 
the  heat  pipe  described  by  Ponnappan  (1989),  with  the.  following  result. 

FSL  ~  1.82  >  1 

Again,  no  problems  were  encountered  by  Ponnappan  during  the  startup  process.  From 
eqns.  (2.26)  and  (2.29),  it  can  be  seen  that  a  larger  working  fluid  inventory  and  larger 
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evaporation  latent  heat  aid  in  frozen  startup,  while  a  larger  heat  capacity  and  higher 
melting  temperature  have  an  opposite  effect.  The  frozen-startup  limit  provides  a  tool 
which  will  aid  designers  in  avoiding  difficulties  during  the  frozen  startup  period. 
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TRANSIENT 
TWO-DIMENSIONAL 
GAS-LOADED  HEAT  PIPE 
ANALYSIS 


3.1  SUMMARY 


A  two-dimensional,  transient  mathematical  model  that  accounts  for  diffusion  and 
variable  properties  on  the  operation  of  a  heat  pipe  is  presented.  The  major  advantage 
over  previous  models  is  that  this  model  treats  the  noncondensible  gas  as  a  separate 
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entity  which  is  described  by  mass  transport  phenomena.  Also,  the  energy  transport 
through  the  wall  is  coupled  to  the  transient  operation  of  the  heat  pipe  through  the  use 
of  a  conjugate  solution  technique.  The  complete  behavior  of  the  heat  pipe,  along  with 
the  location  and  two-dimensional  shape  of  the  noncondensible  gas  front,  are  modeled 
from  the  initial  continuum  flow,  liquid  state  startup  to  steady-state  conditions.  The 
proposed  model  predicted  the  existing  experimental  data  for  the  operation  of  high- 
temperature  heat  pipes  with  and  without  non  condensible  gases. 

3.2  INTRODUCTION 

Gas-loaded  heat  pipes  have  been  applied  in  many  diverse  fields,  and  are  useful  when 
the  temperature  of  a  device  must  be  held  constant  while  a  variable  heat  load  is 
dissipated.  However,  their  use  is  limited  by  the  lack  of  a  complete  understanding  of 
the  performance  of  these  devices,  particularly  during  transient  operation. 

Gas-loaded  heat  pipes  have  been  previously  modeled  with  several  levels  of 
approximation.  The  classical  flat-front  analytical  model  of  Marcus  and  Fleishman 
(1970)  neglected  all  diffusion  across  the  vapor-gas  interface.  Later  studies  modeled 
one-dimensional  steady  diffusion  (Edwards  and  Marcus,  1972),  and  one-dimensional 
transient  diffusion  (Shukla,  1981).  Rohani  and  Tien  (1973)  studied  the  steady-state 
two-dimensional  diffusion  process  in  a  gas-loaded  heat  pipe.  The  importance  of  mod¬ 
eling  radial  diffusion  in  a  gas-loaded  heat  pipe,  was  demonstrated  where  the  noncon- 
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densible  gas  tends  to  accumulate  at  the  liquid-vapor  interface,  which  retards  vapor 
condensation.  However,  Rohani  and  Tien  (1973)  neglected  the  effect  of  conjugate 
heat  transfer  through  the  wall  as  well  as  transient  effects.  The  transient  response 
of  a  gas-loaded  heat  pipe  was  shown  to  be  important  by  Shukla  (1981)  as  transient 
thermal  overshoots  were  discovered.  However,  this  study  did  not  include  the  effects 
of  radial  heat  conduction  in  the  wall  or  wick,  nor  the  radial  diffusion  effects  of  the 
noncondensible  gas. 

Cao  and  Faghri  (1990)  presented  a  two-dimensional  transient  analysis  of 
heat  pipe  dynamics.  The  analysis  included  the  compressible  vapor  flow  and  the 
coupling  between  the  heat  pipe  wall/wick  and  the  vapor  flow.  However,  the  analysis 
was  limited  to  the  conventional  heat  pipe  without  noncondensible  gas. 

In  the  present  study,  the  two-dimensional,  transient  gas-loaded  heat  pipe  op¬ 
eration,  including  the  effects  of  conjugate  heat  transfer  through  the  wall,  is  modeled 
through  a  solution  of  the  general  differential  conservation  equations.  This  proce¬ 
dure  was  used  to  simulate  the  high  temperature  heat  pipe  experimentally  studied  by 
Ponnappan  (1989)  both  with  and  without  noncondensible  gases. 
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3.3  MATHEMATICAL  MODELING 
3.3.1  Vapor  Space 

The  physical  configuration  and  coordinate  system  of  the  gas-loaded  heat  pipe  studied 
is  shown  in  Fig.  3.1.  Gas-loaded  heat  pipes  offer  isothermal  operation  for  varying 
heat  loads  by  changing  the  overall  thermal  resistance  of  the  heat  pipe.  As  the  heat 
load  increases,  the  vapor  temperature  and  total  pressure  increase  in  the  heat  pipe. 
This  increase  in  total  pressure  compresses  the  noncondensible  gas  in  the  condenser, 
increasing  the  surface  area  available  for  heat  transfer,  which  maintains  a  nearly  con¬ 
stant  heat  flux  and  temperature. 

The  differential  conservation  equations  for  transient,  compressible,  two- 
species  flow  with  constant  viscosity  are  as  follows  (Ganic  et  al.,  19S5;  Bird  et  ah, 
1960) 


Mass: 


dp  i  d  ,  d  . 

_  +  +  -ipw)  =  0 


(3.1) 


Momentum: 


-  VP  +  ijiv(  v  •  V)  +  ?v!(? 


(3.2) 
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Figure  3.1:  Gas-loaded  heat  pipe  configuration  and  coordinate 


system. 
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Energy. 


P^~  -  v  •  *vr  -  v  -  ^  Du^TVpj  ~  ~  -  0  (3.3) 

where  subscript  j  denotes  either  vapor  or  gas. 

Species: 


£pg_ 

Dt 


-  V  •  DgvVpg  =  0 


(3.4) 


where 


$  =2 


SHrMt: 


/  dt>  du/ 

+  (&+ 6f, 


1  3  dw 

-—(rv)+  — - 
r  or  oz 


(3.5) 


where  v  and  w  are  the  radial  and  axial  vapor  velocities,  p  is  the  total  mixture  pressure, 
J[  is  the  mass-fraction-weighted  mixture  viscosity,  c p  is  the  specific  heat  of  the  mixture, 
k  is  the  thermal  conductivity  of  the  mixture,  D,i  is  the  self-diffusion  coefficient  for  both 
the  vapor  and  gas  species,  Dqv  is  the  mass  diffusion  coefficient  of  the  vapor-gas  pair, 
pg  is  the  density  uf  the  noncondensible  gas,  and  the  mixture  density  is  p  =  pg  +  p v. 
The  partial  gas  density  is  determined  from  the  species  equation  and  the  vapor  density 
is  found  from  the  ideal  gas  relation  using  the  partial  vapor  pressure. 

The  two  choices  in  species  conservation  formulation  are  mass  and  molar 
fraction.  Molar  fi  .-.tion  offers  the  possibility  of  a  direct  simplification  in  the  formu¬ 
lation  by  the  assumption  of  constant  molar  density.  This  assumption  is  valid  over  a 
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wider  range  of  temperature  and  pressure  than  the  corresponding  assumption  of  con¬ 
stant  mass  density.  However,  when  molar  fractions  are  used,  the  momentum  equation 
must  be  written  in  terms  of  molar- weighted  velocities.  The  resulting  equation  cannot 
be  written  in  terms  of  the  total  material  derivative,  and  is  significantly  more  difficult 
to  solve. 

A  benefit  of  the  general  equation  formulation  is  the  allowance  for  variable 
properties.  Typical  of  compressible  gas  applications,  the  density  is  related  to  the 
temperature  and  pressure  through  the  equation  of  state 


V  = 


•pRT 

W 


(3.6) 


where  M  is  the  molecular  weight  of  the  vapor-gas  mixture,  and  R  is  the  universal 
gas  constant. 

In  the  species  equation,  Dgv  is  a  function  of  pressure  and  temperature.  For  a 
vapor-gas  mixture  of  sodium-argon,  the  relationship  for  Dgl,  is  (Edwards  et  al.,  1979) 


D9V  =  1.3265  x  10_3T3/2(p)~l  (3.7) 

where  T  is  in  Kelvin,  p  is  in  N/m2,  and  Dgv  is  in  m2/s.  Following  a  similar  procedure, 
the  variable  diffusion  coefficient  formulation  for  the  sodium-helium  pair  is  (Edwards 
et  ah,  1979) 
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r 


Dgv  =  1.2795  x  10 -3T3/2(p)_1  (3.8) 

The  interspecies  heat  transfer  whicl)  occurs  through  the  vapor-gas  mass 
diffusion  was  modeled  with  a  self-diffusion  model,  as  described  by  Bird  et  al.  (1960). 
In  the  present  model,  however,  the  self-diffusion  coefficient,  Dd ,  was  assumed  constant 
at  the  initial  temperature  of  the  heat  pipe. 

3.3.2  Wall  and  Wick 

i 

In  the  heat  pipe  wall,  heat  transfer  is  described  by  the  transient  two-dimensional 
conduction  equation 


dT  _  (  dT\  &t_ 

PwCpw  dt  ~  w[r  Or  \V  dr)  +  8z\ 
where  the  subscript  w  denotes  the  heat  pipe  wall  material.  In  the  heat  pipe  wick, 
the  model  uses  the  conduction  simplification  discussed  by  Cao  and  Faghri  (1990). 
This  assumption  states  that,  for  thin  wicks  and  working  fluids  with  high  thermal 
conductivities,  the  effect  of  convection  due  to  the  liquid  flow  in  the  porous  wick  is 


(3.9) 
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negligible.  Therefore,  heat  transfer  in  the  wick  is  by  conduction  only,  as  described 
by: 


dT 

=  keff 


1 md_(.dT\  &T 

r  dr  V  dr  )  +  dz 2 


(3.10) 


where  ( pCp)eff  =  <p(pCp)e  4-  (1  —  <p)(pCp)s,  ip  is  the  wick  porosity  and  ( pcp)s  is  the  heat 
capacity  of  the  solid  wick  structure.  ke ft  is  the  effective  thermal  conductivity  of  the 
wrapped  screen  wick  structure  (Chi,  1976): 


_  kt[(kj  +  ks)  -  (1  -  p){kt  -  ks)} 

'eff  [(A* +  &4) +  (!-<?)(*< -Ml 


(3.11) 


3.3.3  Boundary  Conditions 

At  the  end  caps  of  the  heat  pipe,  the  no-slip  condition  for  velocity,  the  adiabatic 
condition  for  temperature,  and  the  overall  gas  conservation  conditions  are  imposed. 


At  z  =  0 


9T  dPa  n 
V  =  W  =  —  =  =  0 

dz  dz 


(3.12) 


At  2  =  1:  v  =  w=  —  =  0,  Pg  =  pg,BC  (3.13) 

where  p9ibc  is  iteratively  adjusted  to  satisfy  overall  conservation  of  noncondensible 
gas.  This  boundary  condition  is  implemented  through  the  calculation  of  the  total 
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:  lass  of  noncondensible  gas  in  the  heat  pipe.  If  the  total  mass  is  found  to  be  less  than 
the  mass  initially  present  in  the  pipe,  the  boundary  value  is  increased  by  10%  of  the 
previous  value.  Conversely,  if  the  calculated  mass  is  larger  than  that  initially  present, 
the  boundary  value  is  decreased.  This  ensures  the  conservation  of  the  overall  mass 
to  within  a  preset  tolerance,  which  is  1%  in  the  present  formulation.  The  symmetry 
of  the  cylindrical  heat  pipe  requires  that  the  radial  vapor  velocity  and  the  gradients 
of  the  axial  vapor  velocity,  temperature,  and  gas  density  are  zero  at  the  centerline 


(r  =  0): 


dw  dT  dpa  f. 
"=  =  ^=1?  =  ° 


(3.14) 


The  liquid-vapor  interface  (7-  =  Rv)  is  impermeable  to  the  noncondensible  gas. 


mg  —  SsADgyV  Pg  +  PgSgV  =  0 


(3.15) 


where  Sg  is  the  surface  area  of  the  heat  pipe  core.  This  formulation  of  rrig  accounts 
for  both  convective  and  diffusive  non  condensible  gas  mass  fluxes  at  the  vapor-liquid 
interface.  To  ensure  saturation  conditions  in  the  evaporator  section  (and  part  of 
the  adiabatic  section  since  the  exact  transition  point  is  determined  iteratively),  the 
Clausius- Clapeyron  equation  is  used  to  determine  the  interface  temperature  as  a  func¬ 
tion  of  pressure.  The  interface  radial  velocity  is  then  found  through  the  evaporation 
rate  required  to  satisfy  heat  transfer  requirements.  The  no-slip  condition  is  still  in 
effect  for  the  axial  velocity  component. 
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At  r  =  Rv  for  z  <  Le  +  La: 


Tsat  (t0  p o) 

(3.16) 

Vs~  (  \— 

(3.17) 

{hfg  ■+■  CpgT sat  J  pg 

w  =  0 

(3.18) 

where  kg,  Cpg  and  ~ps  are  the  vapor-gas  mixture  properties  at  the  liquid- vapor  interface. 
In  eqn.  (16),  the  saturation  temperature  of  the  vapor  is  found  from  the  partial  vapor 
pressure.  A  solution  of  the  momentum  equation  gives  the  total  mixture  pressure,  but 
the  partial  vapor  pressure  can  be  found  from  using  the  local  gas  density. 

P’  =  KP{'-?)  (319) 

This  equation  was  derived  assuming  a  mixture  of  ideal  gases  following  Dalton’s  model 
for  mixtures. 

At  the  liquid-vapor  interface  in  the  active  portions  of  the  condenser  section, 
vapor  condenses  and  releases  its  latent  heat  energy.  This  process  is  simulated  by 
applying  a  heat  source  at  the  interface  grids  in  the  condenser  section.  The  inter¬ 
face  velocity  can  be  obtained  through  a  mass  balance  between  the  evaporator  and 
condenser  section  allowing  for  inactive  sections  ol  idenser. 

At  r  =  Rv  for  z  >  Lc  -f  La : 

(ho  =  -  {h/g  +  c^Tg)  {p  -  pg)vs  (3.20) 
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Due  to  the  conjugate  nature  of  the  solution  procedure,  the  boundary  condi¬ 


tion  between  the  wick  and  the  wall  is  automatically  satisfied.  This  condition  requires 
the  equality  of  the  heat  fluxes  into  and  out  of  the  wick-wall  interface  (r  =  Rw). 

^f  =  *.wr§  (3.21) 

At  the  outer  pipe  wall  surface,  the  boundary  conditions  depend  on  both  the 
axial  position  and  the  mechanism  of  heat  transfer  being  studied.  In  the  evaporator, 
a  constant  heat  flux  is  specified.  In  the  condenser,  a  radiative  boundary  condition  is 
imposed. 

Evaporator: 


Adiabatic: 


IL! 


dT_ 

dr 


dT 

dr 


Condenser: 


(3.22) 


(3.23) 


(3.24) 
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3.3.4  Initial  Conditions 

Since  this  model  simulates  the  transient  operation  of  a  gas-loaded  heat  pipe,  the  initial 
physical  conditions  of  the  heat  pipe  are  required  for  a  solution  of  the  problem.  The 
initial  temperature  is  uniform  throughout  the  heat  pipe  at  the  start  of  the  simulation. 
There  is  no  motion  of  either  the  gas  or  vapor,  and  the  noncondensible  gas  is  evenly 
distributed  throughout  the  vapor  space  by  diffusion.  The  initial  temperature  of  the 
heat  pipe  is  above  the  free-molecular/ continuum-flow  transition  temperature  for  the 
specific  heat  pipe  vapor  diameter  (Cao  and  Faghri,  1993). 

3.4  NUMERICAL  METHODOLOGY 

The  numerical  procedure  used  in  this  model  is  the  SIMPLE  method  (Patankar,  1980). 
Discretizing  the  differential  equations  resulted  in  a  system  of  algebraic  equations  in 
terms  of  T,  w,  v,  and  pg.  This  system  of  equations  was  solved  using  the  Thomas  or 
tri-diagonal  matrix  algorithm. 

The  numerical  procedure  is  as  follows: 

1.  Input  the  initial  conditions. 

9  Initialize  the  pressure,  temperature,  gas  density,  and  velocity  fields. 

3.  Solve  the  momentum  equations  using  the  mixture  density  to  obtain  v*  and 
w“ ,  which  are  the  vapor  velocities  based  on  the  guessed  pressure. 
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4.  Solve  for  the  new  pressure  field  and  update  the  previous  pressure  field. 


5.  Calculate  v  and  to  from  their  starred  values  using  velocity  correction  for¬ 
mulas. 

6.  Solve  the  energy  equation. 

7.  Solve  the  mass  transport  equation. 

8.  Using  the  mixture  molecular  weight,  determine  the  new  p  from  the  equa¬ 
tion  of  state. 

9.  Steps  3-8  are  repeated  until  convergence  is  reached. 

There  are  two  possible  methods  to  account  for  vapor  compressibility.  The 
first  method  is  to  use  a  density  correction  formulation 

9  =  p"  +  p'  =  Pm  +  A>'  (3.25) 

where  Kp  =  dp/ dp,  and  p'  is  the  guessed  pressure.  The  second  method  uses  the  ideal 
gas  relationship  in  conjunction  with  an  additional  term  in  the  energy  equation.  In 
this  method,  the  pressure  is  chosen  as  the  dependent  variable  and  the  vapor  density 
is  determined  directly.  The  additional  term  in  the  energy  equation  accounts  for  the 
temperature  variation  caused  by  the  high  vapor  velocities,  and  must  be  included  when 
using  the  equation  of  state.  In  the  present  analysis,  the  ideal  gas  relationship  resulted 
in  better  convergence  characteristics.  When  the  vapor  flow  was  subsonic,  it  was  found 
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that  both  methods  gave  very  similar  results. 

In  the  present  study,  a  uniform  axial  grid  and  a  nonuniform  radial  grid  were 
used.  The  numerical  simulation  used  12  (radial)  x  55  (axial)  cells  for  the  2.0  in  long 
heat  pipe.  The  numerical  model  was  run  for  several  different  grid  distributions  with 
less  than  a  3%  maximum  deviation  between  the  solutions.  The  number  of  cells  used 
were  chosen  to  optimize  the  accuracy  of  the  solution  versus  the  required  computer 
time.  The  physical  properties  at  the  interface  are  based  on  the  harmonic  average  of 
the  properties  at  the  adjacent  cells  (Patankar,  1980). 

3.5  RESULTS  AND  DISCUSSION 

Validation  of  the  present  numerical  heat  pipe  model  was  accomplished  by  simulating 
four  cases  of  experimental  data  by  Ponnappan  (1989),  which  are  shown  in  Table  3.1. 
The  heat  pipe  geometry  and  operating  conditions  are  shown  in  Table  3.2.  Vector  and 
contour  plots  for  Case  1  are  shown  in  Fig.  3.2,  where  the  effect  of  vapor  compressibil¬ 
ity  can  be  seen.  Due  to  the  small  vapor  space  aud  high  heat  input,  the  vapor  velocity 
is  very  high  (~  330  m/s).  For  sodium,  this  velocity  corresponds  to  a  Mach  number 
of  approximately  0.5.  From  Fig.  3.2(a),  a  temperature  decrease  can  be  seen  at  the 
exit  of  the  evaporator  section,  where  the  maximum  velocity  occurs.  Furthermore,  as 
the  vapor  velocity  decreases,  a  temperature  recovery  can  be  seen  in  the  condenser 
section.  Also  shown  in  Fig  3.2(a)  are  the  temperature  drop  across  the  wall  and  wick. 
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Table  3.1:  Operating  conditions  for  various  cases 


Case 

Gas  Pressure 

Mass  of  Gas 

Heat  Load 

Condenser 

Number 

(torr) 

(kg) 

(W) 

Emissivity 

1 

0.0 

0.0 

515 

0.3 

2 

1.35 

3.66  xl0~7 

451 

0.56 

3 

1.35 

3.66  xlO-7 

258  to  306 

0.3 

4 

1.35 

3.66xl0~7 

258  to  451 

0.3 
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3.2:  Heat  pipe  geometry  and  parameters 


z  (m) 


0.0000  0.0111 


r  (m) 


0.00  6.35  0.00  6.35 

r  (m)  *10"3  r  (m)  *10"’ 


Figure  3.2:  Numerical  simulation  of  the  high-temperature  heat  pipe  (Case  1):  (a) 
Vapor  temperature  contours;  (6)  Vapor  velocity  vectors  in  the  r-z  plane;  (c)  Vapor 
pressure  contours. 
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The  temperature  profiles  in  the  wall  are  nearly  linear  since  the  thermal  conductivity 
of  the  stainless-steel  pipe  wall  is  relatively  low.  However,  the  direction  of  the  temper¬ 
ature  gradients  in  the  respective  sections  can  be  seen.  As  expected,  the  temperature 
gradient  is  positive  in  the  evaporator  as  heat  is  being  added  to  the  pipe.  Likewise,  the 
temperature  gradient  is  zero  in  the  adiabatic  section  and  negative  in  the  condenser 
section. 

The  vapor  velocity  vectors  in  the  r-z  plane  are  shown  in  Fig.  3.2(6).  Since 
the  evaporator  section  is  much  shorter  than  the  condenser  section,  the  radial  vapor 
velocity  in  the  evaporator  is  much  higher  than  that  in  the  condenser.  The  highest 
vapor  velocity  occurs  at  the  centerline  of  the  heat  pipe,  as  expected. 

The  vapor  pressure  contours  are  shown  in  Fig.  3.2(c),  where  the  two- 
dimensional  effects  are  seen  to  be  negligible.  This  was  expected  since  the  vapor  space 
diameter  for  this  heat  pipe  is  very  small.  More  importantly,  the  vapor  temperature 
and  pressure  are  not  dependent  in  the  vapor  space.  In  a  heat  pipe,  phase  change 
occurs  only  at  the  vapor-liquid  interface,  which  implies  that  saturation  conditions 
are  seen  only  at  the  interface.  In  the  vapor  space,  however,  the  saturation  relation 
does  not  necessarily  apply,  where  the  vapor  temperatures  are  determined  from  en¬ 
ergy  requirements,  accounting  for  viscous  and  compressible  effects.  The  transient 
axial  temperature  profiles  for  the  no-gas  case  are  shown  in  Fig.  3.3  for  Case  1.  From 
Fig.  3.3(a),  the  comparison  with  the  experimental  data  is  very  good.  The  transient 
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1  1.5 

z  (m) 


Figure  3.3:  Temperature  profiles  for  conventional  heat  pipe  with  Qin  =  515  W  (Case 
1):  (a)  Transient  wall  temperature  profile;  ( b )  Transient  centerline  temperature  pro- 


time  for  this  heat  pipe  is  approximately  '2000  s.  The  vapor  dynamics  for  this  case  are 
shown  in  Fig.  3.4. 

The  gas-loaded  heat  pipe  experimentally  studied  by  Ponnappan  was  sim¬ 
ulated  (Case  2),  with  results  shown  in  Figs.  3.5  and  3.6.  This  case  has  a  higher 
radiative  emissivity  than  Case  1,  which  results  in  a  decrease  in  the  thermal  resis¬ 
tance  at  the  outer  condenser  surface.  In  the  experiment  performed  by  Ponnappan, 
the  emissivity  was  increased  when  the  non  condensible  gas  was  added,  so  that  a  near¬ 
constant  operating  temperature  (compared  to  Case  1)  could  be  maintained.  In  Fig. 
3.5,  the  wall  and  vapor  temperatures  decreased  significantly  in  the  con  denser  section 
due  to  the  presence  of  the  noncondensible  gas.  The  gas  density  increased  in  the  con¬ 
denser  during  the  transient  operation,  as  shown  in  Fig.  3.6(a).  The  steady  state  wall 
temperature  is  in  good  agreement  with  the  data  by  Ponnappan. 

To  further  validate  the  numerical  procedure,  the  gas-loaded  heat  pipe  of 
Ponnappan  with  a  low  radiative  emissivity  (Case  3)  was  simulated.  The  velocity 
vector  and  contour  plots  for  the  gas-loaded  case  with  a  low  radiative  emissivity  (Case 
3)  are  shown  in  Fig.  3.7.  In  Fig.  3.7(a),  the  temperature  distribution  for  the  gas- 
loaded  case  is  significantly  different  from  that  of  the  no-gas  case.  As  expected,  there 
is  a  large  temperature  drop  across  the  inactive  condenser  section,  and  the  radial 
temperature  variations  in  this  section  are  relatively  uniform.  This  is  due  to  the  low 
vapor  velocities  in  this  region  and  the  resulting  low  convective  heat  transfer  in  the 
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900 


Figure  3.5:  Temperature  profiles  for  gas- loaded  heat  pipe  with  Qin  =  451  W  (Case  2) 
(a)  Transient  wall  temperature  profile;  (6)  Transient  centerline  temperature  profile. 
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Figure  3.7:  Numerical  simulation  of  the  gas- loaded  high-temperature  heat  pipe  (Case 
3):  (a)  Vapor  temperature  contours;  ( b )  Noncondensible  gas  density  profiles;  (r) 
Vapor  velocity  vectors  in  the  r-z  plane;  (d)  Vapor  pressure  contours. 
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inactive  condenser  section. 


The  vapor  density  contours,  shown  in  Fig.  3.7(6),  are  radially  uniform 
and  have  a  sharp  axial  gradient.  This  is  a  conduction-controlled  heat  pipe,  with 
axial  conduction  through  the  heat  pipe  wall  and  wick  as  the  principal  heat  transfer 
mechanism  across  the  vapor-gas  interface. 

The  results  for  Case  3  with  a  heat  input  of  258  W  are  shown  in  Figs.  3.8 
and  3.9.  In  Fig.  3.8(a),  the  steady-state  wall  temperature  is  compared  with  the  ex¬ 
perimental  data  by  Ponnappan  (1989)  with  reasonable  agreement.  In  Fig.  3.9(a),  the 
transient  gas  density  profiles  are  given,  where  the  importance  of  including  conjugate 
heat  transfer  with  the  wall  can  be  recognized.  If  the  gas  distribution  were  controlled 
purely  by  the  vapor  dynamics,  the  transient  time  would  be  on  the  order  of  10-3 
s  (Shukla,  1981).  Due  to  the  presence  and  subsequent  damping  effect  of  the  wall, 
the  transient  time  of  the.  noncondensible  gas  distribution  is  substantially  increased. 
Physically,  this  phenomenon  is  controlled  by  temperature-dependent  diffusion.  At 
the  beginning  of  operation,  after  liquid-  or  solid-state  startup,  the  vapor  temperature 
in  the  condenser  is  high.  This  high  temperature  increases  the  diffusive  mass  flux 
across  the  vapor-1  as  interface.  However,  as  time  increases,  the  vapor  temperature 
in  the  condenser  decreases,  thus  decreasing  the  mass  flux  of  the  noncondensible  gas 
across  the  interface,  and  increasing  the  partial  gas  density  in  the  reservoir.  This 
phenomenon  is  reinforced  by  the  compressible  behavior  of  the  vapor.  As  shown  in 


73 


O  t  =  50.0  s 

A  ♦  =  100.0  s _ 

+  t  a  500.0  8 _ 

x  t  =  1000.0  s 
O  t  =  2000.0  » 

V  t  a  3000.0  r 
8  St»ody-Slot« 

•  Ponnappon  (1989) 


V^K 


Figure  3.8:  Temperature  profiles  for  the  gas- loaded  heat  pipe  with  Qin  =  258  W 
(Case  3):  (a)  Transient  wall  temperature  profile;  ( b )  Transient  centerline  temperature 


Fig.  3.9(5),  the  vapor  velocity  is  initially  low.  Therefore,  the  inertial  effects  of  the 
vapor  against  the  stationary  noncondensible  gas  are  small.  As  time  increases,  the 
vapor  temperature  increases  in  the  evaporator,  causing  the  vapor  density  to  decrease. 
Thus,  to  maintain  the  same  mass  flow  rate,  the  vapor  velocity  must  increase.  This 
velocity  increase  causes  an  increase  in  the  inertial  forces  acting  on  the  noncondensible 
gas,  which  further  compresses  the  gas  slug. 

After  the  gas-loaded  heat  pipe  reached  the  steady  state  for  a  heat  input 
of  258  W,  the  heat  input  was  pulsed  to  306  W,  as  shown  in  Figs.  3.10  and  3.11. 
The  transient  axial  temperature  profiles  a-e  shown  in  Fig.  3.10,  and  the  transient 
vapor-gas  dynamics  are  shown  in  Fig.  3.11.  The  curves  labeled  with  t  =  0.0  s.  are 
corresponding  to  the  steady  state  for  Q,„  =  258  W.  Of  primary  concern  in  Fig.  3.10 
is  the  comparison  of  the  transient  time  for  the  258  to  306  W  heat  pulse  with  that  of 
the  initial  response  to  the  258  W  heat  input.  As  expected,  the  pulsed  transient  time 
is  nearly  30%  of  the  initial  response  time.  Physically,  this  is  due  to  the  conduction- 
controlled  behavior  of  the  heat  pipe.  In  the  initial  258  W  response,  the  heat  stored  in 
the  condenser  pipe  wall  needed  to  be  rejected  to  the  ambient  as  the  noncondensible 
gas  plug  developed.  However,  at  the  application  of  the  258  to  306  W  heat  pulse,  the 
condenser  temperature  was  much  closer  to  its  final  value,  which  required  much  less 
time  to  adjust. 

As  the  operating  power  changed  from  258  to  306  W,  the  vapor  pressure 
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Figure  3.10:  Temperature  profiles  for  the  gas- loaded  heat  pipe  with  a  pulsed  heat 
input  of  Qin  =  258  to  306  W  (Case  3):  (a)  Transient  wall  temperature  profile;  ( b ) 


Transient  centerline  temperature  profile. 
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Figure  3.11:  Vapor-gas  dynamics  for  the  gas-loaded  heat  pipe  with  a  pulsed  heat 
input  of  Qin  =  258  to  306  W  (Case  3):  (a)  Transient  centerline  gas  density  profile; 
(6)  Transient  centerline  axial  velocity  profile. 


of  the  heat  pipe  increased.  This  pressure  increase  compressed  the  nonc.ondensible 
gas  plug,  and  increased  the  active  condenser  length,  as  can  be  seen  in  Fig.  3.11(a). 
However,  this  pressure  increase  caused  the  vapor  density  to  increase,  due  to  the 
vapor  compressible  effects.  As  the  density  increased,  the  vapor  velocity  decreased  to 
maintain  a  constant  mass  flow  rate.  As  before,  these  compressibility  effects  can  oe 
seen  in  Fig.  3.1 1  (Z>). 

The  transient  axial  temperature  results  for  a  pulsed  heat  input  of  258  to 
451  W  (Case  4)  are  shown  in  Fig.  3.12,  and  the  vapor-gas  dynamics  are  shown  in 
Fig.  3.13.  The  main  factor  to  be  examined  from  this  case  is  the  initial  response  of 
the  gas-loaded  heat  pipe,  as  shown  in  Figs.  3.12(6)  and  3.13(6).  As  the  heat  input 
is  pulsed  from  258  to  451  W,  the  vapor  temperature  initially  dropped  at  the  inlet 
to  the  adiabatic  section.  Similarly,  the  axial  velocity  at  that  point  increased  nearly 
60%.  These  effects  have  been  observed  by  Cao  and  Faghri  (1990)  and  are  attributed 
to  compressibility  effects. 

3.6  CONCLUSIONS 

A  transient  two-dimensional  simulation  of  gas-loaded  heat  pipes  has  been  performed. 
The  present  numerical  model  used  a  conjugate  solution  technique  to  couple  the  un¬ 
steady  heat  transfer  in  the  wall  and  wick  with  the  transient  vapor  How.  The  boundary 
condition  at  the  outer  condenser  pipe  wall  influenced  the  operating  temperature  of 
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Figure  3.12:  Temperature  profiles  for  the  gas-loaded  heat  pipe  with  a  pulsed  heat 
input  of  Qin  =  258  to  451  W  (Case  4):  (a)  Transient  wall  temperature  profile;  ( b ) 
Transient  centerline  temperature  profile. 
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Figure  3.13:  Vapor-gas  dynamics  for  the  gas-loaded  heat  pipe  with  a  pulsed  heat 

input  of  <2i„  =  258  to  451  W  (Case  4):  („)  Transient  centerline  gas  density  profile; 
(b)  Transient  centerline  axial  velocity  profile. 
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the  heat  pipe  and  altered  the  position  of  the  vapor-gas  interlace.  Increasing  the  ra¬ 
diant  emissivity  decreased  the  operating  temperature  and  transient  time  period  of 
the  heat  pipe.  Steady-state  gas-loaded  heat  pipe  operation  depends  heavily  on  the 
mass  of  noncondensible  gas,  as  this  parameter  affects  the  operating  temperature  and 
active  condenser  length.  It  was  also  shown  that  the  transient  time  period  increased 
when  gas  was  added  because  more  heat  needed  be  conducted  through  the  pipe  wall 
to  the  previously  active  part  of  the  condenser  section.  While  the  noncondensible  gas 
is  fundamentally  one-dimensional  in  nature,  the  complete  problem  must  be  solved  in 
two  dimensions  accounting  for  conjugate  heat  transfer  through  the  wall  and  wick. 
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